APPROXIMATION OF SCHRODINGER UNITARY GROUPS OF 
OPERATORS BY PARTICULAR PROJECTION METHODS 



FREDY VIDES 



Abstract. In this paper we work with the approximation of unitary groups of 
operators of the form e~ ttH where H £ ^(H) is the self-adjoint Hamiltonian 
of a given Hermitian quantum dynamical system modeled in the discretizable 
Hilbcrt space H = 'H(G), to perform such approximations we implement some 
techniques from operator theory that we name particular projection methods 
by compatibility with quantum theory conventions. Once particular represen- 
tations are denned we study the interelation between some of them properties 
with the original operators that they mimic. In the end some estimates for 
numerical implementation are presented to verify the theoretical discussion. 



1. Introduction 

In this work we will focus our attention in the approximation of Schrodinger 
semigroups of operators that will in be described in general by the set {Ut := e~ ltH : 
t € R}, whose elements clearly satisfy the semigroups conditions: (i) Uq = 1, (ii) 
Ut°U s (-) = U t + S {-) and (hi) lim /l _ ! . + UhX — x, Vx G V C H(G), besides by theorem 
T IA.3l the condtion (iv) ||C/(Mo|| = ||ito|| ,i G R, will be also satisfied when H is self 
adjoint, wich means that Ut is unitary for any t £ R. In the expression presented 
above the operator H 6 J>f (%) is the Hamiltonian of a given quantum dynamical 
system whose abstract evolution equation is given by: 

Eip(t) = Hi>(t) (1.1) 

with ^(0) = -0o £ H and where E — > j^dt, here H g Jzf(%) will in general have 
the form H — p^p + V(-) with pt — > _j_v + b and with V E C a=1 {B C H), for 
simplicity, in this work we will consider our scale such that % = 1, also we will 
have that in some suitable sense the operator H £ ^(H) will be restricted by some 
boundary conditions related to the media where a particular quantum dynamical 
system evolves. 

In the following sections we will implement some operator theory techniques in 
the theoretical analysis of the approximation schemes of the Schrodinger semigroups 
and in the end some numerical implementations will be presented. 
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2. Basics of Quantum Dynamical Systems 

Quantum dynamical systems can be studied using several types of operators, in 
this work we will consider in general quantum dynamical systems that evolve in a 
particular discretizable space of states 71(G), i.e., a separable reproducing kernel 
Hilbert Space with G CC R . A vector tp(t) € H(G) that satisfies (jl.ip receives 
the name of wave function, the evolution of a given quatum system can be described 
using the wave function, the time evolution of the wave function can be computed 
using the corresponding Schrodinger semigroup related to a particular quantum 
dynamical system using the following expression 

V>(i) = e -*"Vo. (2.1) 

The vector tp(t) € H(G) can also be considered like a probability amplitud related 
to measurements concerning to the position of the particles in a quantum system 
modeled by an abstract Shrodinger evolution equation like Usually we can 

express quantum evolution equations using the Dirac's braket notation, we usually 
have that the ket operation |-} : H — > Ti, is defined explicitly by 

|0):=0>eH, (2.2) 

on the other hand we have that the bra operation (-| [■] : H x % — > C, is defined by 

{<f>\[-]:=(;4>),<p€U, (2.3) 

it can be seen that the natural inner product of the space of states T-L can be 
expressed using the Dirac's braket notation in the form 

for any quantum operator A G Jzf('H) and any pair £,0 6%, the operation (A£, <fr) 
can be expressed explicitly by 

(M, 4>) = (01 [M] = (01 a |0 = (4>\M) , e, g n. (2.5) 

The probability density for a specific time will be given by IV 7 ^)! 2 = V'WV'Wj * ms 
statistical approach of the wave function, and the corresponding inner product of 
H(G) alows us to compute the expectation (B) t S R of a given observable B, i.e. 
a quantum operator B G using the following expression: 

w *- mim)u {G) [ ' 

wich in Dirac's braket notation is equivalent to 

= mm) ) • (2 - 7) 

For any two operators A, B 6 «Sf('H(G)) on discretizable quantum space of states 
11(G) we can define a commutator operation by 

[A, S] := AS - (2.8) 

we will have that two operators 1,7 6 _§f (H(G)) commute if, and only if [A", Y] = 
0, when two operators commute we also say that they are compatible. In general 
the operation [■, ■] : J*?(H) x JC(H) — > H defines an operator in jSf (U), If we obtain 
a nomalized representation |\&(t)) := \i>(t)) / \\ip(t)\\ of the wave function ip(t) E H, 
then (|2.7p can be represented by 

(B) t = (*(t)\BMt))- (2-9) 
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Theorem 2.1. Constants of Motion. If a quantum operator A G Jz?(%) is constant 
in time and compatible with the Hamiltonian H 6 J?f ("H) in (jl.ip . and if H G J£{Ti) 
is symmetric then {A) t defines a constant of motion. 

Proof. Since H G ^(H) si symmetric, it can be seen that 

j t {A) t = i(*(t)\HA\*(t))-i(*{t)\AH\*(t)) 

= i(9(t)\[H,A]\*(t)) 

since we also have that A and H are compatible, then [H, A] = and ^ = 0, 
therefore (A) t is a constant of motion. □ 

3. Particular Projection Methods 

In this section we will describe the approximation techniques implemented for 
spatial discretization of the operators related to the dynamics and physical mea- 
surements of the quantum systems described here. 

3.1. Particular Projection in H(G). In this section and some other below we 
will work with the spatial discretization of operators that are present in a quantum 
dynamical system in Schro-dinger picture. When we want to build a discretization 
of a given spatial operator we first need to define a grid, wich is set G m ft C G C R w 
that depends in some suitable sense of the parameters m G Z + and h G R + , in 
particular the cardinality of the grid denoted by N m> h '■= \G m ^h\ depends on m, h 
throug the following rules N m ^h > N m > t h, m > m! and N m j t > N mt h' , h < hi . 

Once we have defined a grid G m ^h C G on the media where a quantum dynamical 
system evolves, we can define a particular projection Pg G J^ > (V.(G),Hs{G)) for 
Us ^ H, with respect to this grid, wich is a projection that can be factored in 
the form P$ := psp^, where ps G J^(H' s ,Hs) and p$ G J£{'H,'H' S ) are particular 
summation and decompostion operators respectively. A particular decomposition 
operator p* s G J2?(H(G),'Hs(G)) is related to a given grid G m ^ C G through the 
following explicit definition 

P f s v := {c k (v,G m , h )} G C Nm - h , v G H(G) (3.1) 

the particular summation operator is related to a given set H(G) D B D B m := bk 
that is a subset of a particular basis B C "H(G), whose elements are compatible 
in some suitable sense with boundary conditions of wave functions in the physical 
media, and that satisfies \B m \ = \G m> h\, thorugh the following expression 

p s x:=Y,Xkbk,x = {x k }&C N ™> h (3.2) 

the functionals Ck(-, G m .h) used for the definition of particular decomposition oper- 
ators are considered in general to satisfy that for any b G B m , psPgb = b, sometimes 
we also impose the condition Ck(bj, G m ,h) = a j^k.j, with 6i j the Kronecker delta, 
the above condition is called pseudo-orthogonality or pseudo-orthonormality when 
ctj = l,Vj. For a given particular projection Ps G J??(H,Hs) we say that it has 
approximation order v mi where v m is a number that depends in some sense on the 
grid paremeter m, if we have that for any v G T~L(G) the projection satisfies the 
relation 

\\P s v - v\\ < c v h^ (3.3) 
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with respect a prescribed norm ||-|| in 7i(G) and where h is the mesh size of the 
prescribed grid. Sometimes we represent a paricular projection Ps £ Jzf (H, Tis) in 
an alternative form given by P m .h where m and h are the grid parameter and mesh 
size respectively. 

3.2. Inner Product Matrices. For a given discretizable Hilbcrt space H(G) 
whose inner product is induced by the inner product map A4 : H x H — > C in 
the form 

(u,v) u :=M[u]{v) (3.4) 

we can obtain a discrete representation M§ of the inner product map M. with 
respect to a particular projection P$ S J??(H,Hs), through the follwing explicit 
definition 

(•)*Ms[.] := M[p s -](p s -) (3.5) 

the matrix Mg receives the name of inner product matrix. The followng result was 
proved in [14] . 

Theorem 3.1. Every inner product matrix is symmetric positive definite. 

If for any u, v £ H we take u := pgU and v := p s v, it can be seen that we can 
express the operation (Psu, Psv) using the inner product matrix in the following 
way 

(P s u, P s v) := M [psP^u] (pspIv) = v*M 5 u (3.6) 

since Wis is symmetric positive definite, and if we denote by Wg the formal square 
root of Ms , we can obtain an alternative expression for the above operation that 
will be given by 

(P s u, P s v) = (W 5 v)*(W s u) = (W 5 u,W 5 v) 2 (3.7) 

where (•, -} 2 is the complex euclidian inner product. In a similar way one can 
express the operation ||P,su||, where ||-|| is the norm induced by (•, •) in H through 
the following relations 

||P 5 t*|| := ((Psu, Psu)) 1 / 2 = ((Wsi^Wsuy 1 / 2 = ||W 5 u|| 2 . (3.8) 

We can express kets and bras in a discrete frame with respect to a particular 
projection Ps £ ££(^K), using the rules 

Wl ^>|*}:=* (3.9) 
(ip\ -^U (*| 4>*M (3.10) 
for ket and bra operations respectively. 

3.3. Particular Representation of operators. We can obtain discrete repre- 
sentations of operators in j£f (%,%') for dicretizablc Hilbcrt spaces, with 
respect to particular projections Ps G Jzf('H) and Qy € Jz?(%'), the correspond- 
ing discretization will receive the name of particular representation, for any given 
A E J^?(H,H') we denote its particular representation with respect to Ps and Qy 
by A 6 C WvXiV ' s and we define it explicitly in the following way 

A := qlAps (3.11) 
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A particular representation A £ C NvxNs is said to have approximation order /j, m , 
with [i m a number that depends on a given grid parameter m related to the par- 
ticular projections, if it satisfies the relation 

\\Q v AP s u-Q v Au\\= W v (Au -q^Au) ^< c u h^ m (3.12) 

with respect to a prescribed norm ||-|| in W , where h is the mesh size of the pre- 
scribed grid, an alternative expression for this property can be obtained when 
% = %' and Ps = Qvi taking ca ■= sup M c u , for c u in the above equation, and then 
writing: 

\\[Ps,A]\\<c A h^ (3.13) 
when h — > + we say that the pair P$, A e almost commute. 

Remark 3.1. It is important to note that for any given B E Jtf(H) and any 
particular projector P$ e J?(H,'Hs), we will have that Ps4> = <f>, V</> G Hs '■= PsT~i 
and B4> = P s B4>,V(f) e U s . 

3.4. Exactly Factorizable Operators. Given two Hilbert spaces X, Y, an oper- 
ator A : X — > X is said to exactly factorizable, if it can be writen in the form 
A := BC, with C : X — > Y and B : Y — > X, in this article we will focus our 
attention on exactly factorizable operators of the form A := BC, with B^ := aC, 
for a = ±1 E R\{0}, these conditions imply that 

(Au, v) x = {Cu, B^v) y = a {Cu, Cv) Y (3.14) 

wich permits us to obtain the following relation 

{Au, v) x — a {Cu, Cv) Y = (u, aC^Cv^ y = (u, BCv) x = {u, Av) x (3.15) 

wich implies that exactly factorizable operators of this type are self-adjoint. Now if 
we take two finite rank particular projections Ps £ ^{X, X£) and Qy G _£f (Y, Y$), 
and if we define the particular representations A := p^Bps, B := p^Bqy and 
C := q^Cps for A, B and C respectively, we can first note that 

V*M V CU = {CU, V) v = (U, CV) 5 = V*(C f )*M s U (3.16) 

wich implies that := M^ 1 C*My, then from the definition of particular pro- 
jections we can obtain the following relations 

V*M S AU = {AU,V) S (3.17) 

= a{CU,CV) v (3.18) 

= a(CV)*M v CU (3.19) 

= V*C*M v [aC]U (3.20) 

The above expressions permit us to represent A by A := M^ 1 C*My[aC] = aC^C, 
and from this we will have that the particular representation of A preserves self- 
adjointness and operator sign (positive/negative) according to a e R\{0}. 
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4. Discrete Time Integration and Matrix Schrodinger Unitary 

Groups 

As we discussed in iJTJ we will consider that the Hamiltonians H G Jz?(H(G)) in 
Scrodinger models like (jl.ip . have the form H = Hq + V(-), where Hq defined by 

H a :=p^p (4.1) 

and with p G &(H,H') and pt e £'(U','H), we will consider that V G C a=1 (V), 
DCfinH'. The abstract initial value problem 

w(t)) = -iH \m)-iv(\m)) U9) 
\m) = m [ ] 

derived from (jl.ip , can be rewrited introducing the integrating factor e ttH ° , in the 
form 

|^(t)) = e - ltHo |^o> + / e^ s - t)Ho V(\iP(t)))ds (4.3) 



it can be seen that Hq is exactly factorizable, hence, self-adjoint and by theorem 
T IA.3I we will have that e ±ltHa ,t G R will be unitary, since we also have that 
V G C Q=1 (£>), we can derive the following result. 

Theorem 4.1. TTie operator U G SSiJ)) defined by 

U(cf>) := e- l * ffo IV'o) + / e^'-^Vfl^t))^ (4.4) 



is a strict contraction with respect to the norm W'Wqh-tt] v(g)) defined by 



U\\C(\-TT] 25(G)) : = SU P U\\v(G) G V ( G )- ( 4 - 5 ) 

te[-T,T] 



when T < 1/cy. 



Proof. Since for any t e R e ± ltH o [ s un itary and since V G C Q_1 (2?) we will have 
that 

11^(10)) - U (\^))\\c([-T,T],T>{G)) = 



e^-')^[U(|0(t)))-U(|a;(t)))]d 



s 

t 



< 

C([-T,T],X>(G)) 



|||0(i)> - k(«)}|lo([-r,ri,p(G)) / ds < 

Jo 

<vT[[|^(i))-| W (t)}|| 

C([-T,T],D(G)) 

taking T < 1/cy the result follows. □ 

Corollary 4.1. There exists a unique solution for (|4.2I) in a time interval [—T,T] 
C R, where T G R + satisfies the restriction presented in the theorem above. 

Proof. Follows from TED □ 
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4.1. Matrix Schrodinger Unitary Groups. If we consider that the space of 
states is discretizable, then we can compute a particular representation Ho of Hq, 
of approximation order v m , that is exactly facotrizable since Hq can be exactly 
factored, hence self adjoint, and will be the generator of a matrix unitary group 
{U := e~ ltH °}. It can be seen that the successive approximation method can 
be performed in two stages in order to compute an approximate discrete solution 
to (|4.2[) in time and space. First we can obtain a generic representation of the 
approximate expression of particular elements of the discrete Schrodinger semigroup 
{Uk ■= e ™ D , h 6 R + , k E Z}, this can be done using a Picard type method in 
the following way. We start with the problem 

i|* , (t))=H |*(t)> (4.6) 

with |^(0)) = l^o), using the succesive approximation method with approximation 
operator defined explicitly by 

T(tf (t)) := |tf„> + / Ho|*(t)> (4.7) 
Jo 

we can obtain the approximatin sequence defined in D IA.2I explicitly given by 

" 1 

\* n (h)) := U |*„> = E M HTH ° )fc '* o) (4 - 8) 

fe=0 

where r 6 R + satisfies the Picard's restriction r < HWnHIoWg 1 1| . It can be seen 
that for these types of initial value problems the succesive approximation method 
produces a solution that coincides with the Taylor expansion of the exponential 
matrix e~ lTM ° , we can take adavantage of this using the Pade approximant of e~ iTM ° 
in the interval [0,r] C R for r € R + a basic time step, denoted by U € Mjv m h (C) 
and defined by U := i? pp (— ztHq),p := [n/2\ where R PP {— irMo) is defined by: 



R pp (-iTU ) := J D pp (-irH )~ 1 iV pp (-zTHo) (4.9) 



with 



N pq (^rm Q ) := E k + * A -irM ) (4.10) 

D p9 (-zrH ) := E /fj-^^l ^) ^ 

f-^ (p + q){r-{i - jV- 

It can be seen that taking § := Npp(— irMo), we will have that D pp (— irHo) = S', 
and if we take § + := (§^) _1 , then we can express (|4.9[) in the form 

U = §+§ (4.12) 

From the relation of (|4.9[) with the Taylor expansion and succesive approximation of 
e -irH and tiie pi car d' s resctriction for r = h T \\MqW~ 1 < HHollg 1 , with < h T < 1 
in (|4.8p . we can obtain the following estimate 



Lemma 4.1. ||e~ lrU ° - U|| 2 < 



(2p+l)! C P,2 P +1 



h 2 P +1 
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Proof. 



-irW 



< 



< 



< 



oo _^ 

E (^- c P ,fe)(-^ H o) 



fe=2p+l 
1 



(2p+l)! 
1 



(2p+l) 



— c p,2p+l 

— c p,2p+l 



Po||^ 



(4.13) 

(4.14) 

(4.15) 
□ 



Now if we want to compute the approximation of e JtH ° corresponding to the 
interval [0,t] C R, with t := mT,m 6 Z+, wc will have that e - itH ° = ( e ^ rH o) m 
becomes approximated by U m , giving this an explicit definition for the discrete 
unitary group {Uk := V h ,k € Z+}. Since H will be considered in general self 
adjoint, i.e., Hj = H , we will have that H is normal, hence can be factored in the 
form H = VH>oV*, with D := diagje^}, and taking A := R pp {— ztO ) we obtain 

3+S 



u = 



= VA+AqV* 



(4.16) 
(4.17) 
(4.18) 



wich implies the following result. 
Lemma 4.2. U*U = 1. 

Proof. 

ITU = 



VASA^VVA^AoV* 

VA^A^A+AqV* 

VASA+Aq^oV* 



from the above relations can see that the operator 

U := Wf u 

satisfies the following relations 



wich implies that 



jW^UWh^L = ||UW H *| 



1*11 



(4.19) 
(4.20) 
(4.21) 
(4.22) 
(4.23) 
□ 

(4.24) 
(4.25) 



= 1, the adjoint of U can be obtained in the following way: 



= Mi 1 (Wi 1 UW H )*M i 
= M H 1 WJjU*W+M H 
= W5 1 W+W^U*W+Wi 



Wi 1 u*w H 



(4.26) 
(4.27) 
(4.28) 
(4.29) 
(4.30) 
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the discrete Schrodinger unitary group relative to Ho, will have the form {Ufc := 
U fc (-),fc G Z + }, wich is consistent with the normalization presented at It can 
be seen that 

UtU = 1 = Ml* (4.31) 
and this implies that the discrete time reversal Schrodinger unitary group will be 
given by {U_fc := (tr) k (-),k £ Z + } and will be coherent with the local time 
reversibility of Schrodinger unitary groups. From lemma L 14.ll we can obtain the 
following. 



Lemma 4.3. 



-imrHo £jr. 



< 



,2p+l 



(2p+l)!"r 



Proof. If we denote by c p ^ the multinomial Pade coefficients, then we will get 



-imrMo 



w 



E 

fc=0 



(-imrHo) fe 



< 



< 



< 



fe! 

X i. 

E, m 

fc=2p+l 
m 2p+l 



k=0 



C p , k )(-lTW ) k 



(2p+l)! 

m 2p+l 



Cp,2p+1 



2p+l 



||H 



i2p+l 



(2p+l)! r 



(4.32) 

(4.33) 

(4.34) 

(4.35) 

□ 



If for a given self-adjoint operator i/i £ J^?(H) we compute its representation 
with respect to a particular projection P$ denoted by Hi, then we can obtain the 
following convergence result. 



Lemma 4.4. 



< Cl h^ + c 2 h 2 T P +1 . 



Proof. By the properties of particular projections we get that 



-irnrHi 



< 



-imrHi 



< c Hl h Vn 



- e 

m 2p+l 



h 2p 

(2p + l)\ T 



+ i 



and taking ci := ch 1 and ci 



(2p+l)< 



the result follows. 



(4.36) 
(4.37) 
□ 



Now if for a particular representation Ho of Hq we denote by {U^ , k £ Z + } and 
by {U_fe,fc G Z + } the direct and reverse Schrodinger semigroups respectively, we 
can obtain a fully discrete representation for (|4.4[) in the form 

k 

:=U fe |*o}+E w ^ T ) 1& ^ y (-) ( 4 - 38 ) 

in this expression {^(r)} are chosen such that the integration in time is exact for 

the polynomial structure of the Taylor similar part of U, wich allows to mimic 
the operator defined in (|4.4j) . Using the properties of particula projection methods 
we can derive the following discrete results. 



10 



FREDY VIDES 



Theorem 4.2. The operator is a strict contraction with respect to the norm 
IHIc([-t,t],Z>(G))> when V e C a=1 {V C H) and h < 1/cy. 

Proof. First it is important to note, that since the quadrature rule defined by 
the sequence {wj(r)} is exact for the polynomial structure of U and if we take 
vp- 7 := ^(jr/k) for any given function ^ e C([~t 1 t],T-L{G)), we will have that 



S fc (*(t))-Sfc(T(t)) 



< 



u fc . 



- v(T J ')| 



< ^^ fc (T)cy || & -TP | 
i=0 



the last relation implies that 
§*($(*)) -S fc (T(f)) 



C([-t,t],XI(G)) 



sup 

te[-T,T-l 



§*(*(*)) -§ fc (T(t)) 



< c vsup||^-T^|| H 
V=o / ' 

< c v t ||$ -T|| c([ _ TiT])X)(G)) 

and taking r < 1/cy the result follows. □ 
Corollary 4.2. T7ie semidiscrete Schrodinger initial value problem 

*i* / (t))=Hoi*(i))+ni*(*)» f4 39) 

|*(0)) = *o 1 ' ; 

/las a unique solution determined by |^q) on cm interval [— t, t] /or r < 1/cy. 

iYoo/. Follows from T|XTJ □ 

Lemma 4.5. ^4 solution \^{t)) to (|4.39[) satisfy the following estimate with respect 
to the its appoximating sequence 

lll*(<)> l*n(*)>[lc([-r,r], C (G)) < j _ ( 2 II I *0> II C ([-r,r], B (G))+^r) (4.40) 

w/iere My := sup,,. ||V"(x)||h. 
Proof. It can be seen that 



IPi(t))-|*o(t)>||h = S(|*o»-|*o> 



Ufc |*o) +£i«i(A!)Ufc_jV(|*o» - |*o) 



3=0 
k 



(4.41) 



(4.42) 



< 2\\\^o)\\ u + J2 w ^ 

< 2|||*o>|| h + t||V(|*o»|| 



u 



k-j 



||^(|*o»|| H (4.43) 
(4.44) 



APPROXIMATION OF SCHRODINGER. UNITARY GROUPS 



11 



also it can be seen that sup T |||*i(r)) - |*o(r))|| H < 2 || |* ) \\c([-t,t],v(G)) + Myr, 
replacing above expressions in R IA.ll and taking r < 1/cy and K — cyr the result 
follows. □ 

Theorem 4.3. If we denote by \ip(t)) £ T-L and by \^ n (t)) £ H the exact solution 
to (|4.2p and the n-th time approximating solution to (|4.39|) respectively and if we 
have that particular representations of quantum operators in H have approximation 
order v m . Then we will have that \^(t)) £ H and |4 , „(i)) £ H satisfy the following 
estimate 

(cvr) n 

l*»(*)>llc([-r,r],7>(G» ^ C ^ m+ l ~ (c V T) {2 ^^C([-r,r h V(G)) + M V T ) 

(4.45) 

Proof. If we denote by |\&(t)) £ H the exact solution to (|4.39j) then we will have that, 
since particular representations of quantum operators in H have approximation 
order v m , we will have that \\\ip(t)) - |*(t))|| H < C^,{t)h Vm ,Vi £ [-r,r],0 < r < oo, 
taking C0 := sup t C^(t), we will obtain the relation || \ip(t)) — \^(t)) \\c([-t t] v(g)) — 
c^h" m . Hence 

\\m)) - l*n(*)>llc{[-T,r],7>(G» ^ H^W> " l*( t ))Hc([-r,r],7>(G)) 

+ \\\*(t))-\* n (t))\\ c([ _ TiT]MG)) 

~ C *^ + l3^( 2 IH*»>llo([-r ( r],l>( C )) 

+M v t). 

wich provides the desired result. □ 

Remark 4.1. Remark R \3.1\ imvlies that for any two compatible operators A 7 B £ 
_£? (H) we will have that 

O = P s 0P s = P S [A,B}P S (4.46) 

= Ps{AB-BA)P s (4.47) 

= P S ABP S - P S BAP S (4.48) 

= P S AP S BP S - P S BP S AP S (4.49) 

- P S AP 2 S BP S - P s BPlAP s (4.50) 

= AsB,s-B,s4s (4.51) 

= (4-52) 

Corollary 4.3. Discrete constants of motion. A particular projection A £ jSf (H) 
o/ an operator A £ that is compatible with Hq £ JC('H) defines a constant of 

motion (A) t for the system modeled by: 

' id t |*(i)> = Hq |*(f)) 



|*(0)> = *o (453) 
where Ho := P^Hqps and Ho := p^Hops. 

Proof. From R0we get that [A,H ] = O, then the result follows from TOJ] □ 

Remark 4.2. 7i is important to observe that [A m ^-, -B m ,/i] = O, V/i G R+, i/ien we 
wiZ/ /iai>e i/iai [A, £?] = 0, with A := lini/ ! _ ! .o+ -^-m,/t a7lc ^ ^ : — hm/ l _ ! .o+ B m ^. 
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5. Examples 

In this section we will present some examples of particular implementation of 
approximation schemes and some related estimates. 

Example 5.1. Quantum Harmonic Oscillator. Given a quantum system on a 
square box G — [0, L] x [0, L], L < oo whose time evolution is described by a quantum 
harmonic oscillator equation that has the form 

l ^- t \m)--=&2\m) + (x 2 +y 2 )\m) (5.1) 

with ip(0) — ipo G <5(G) C H(G)and if the model is restricted to physical condi- 
tions that prevent any representative particle to be located an the boundary, i.e., 
I IV'(^)) | 2 = 0,x G dG, t G R. Physical restrictions imply null Dirichlet boundary 
conditions for (|5.1[) . If we take Hq := A2 := V • V = p^p, wich is clearly exactly 
factorizable, and if we take its particular projection to H a spectral element space 
of first order, then we can express (|5.ip in the form 

l Jt ^ := M ~ lp * M ' p !*(*)) + ( x2 + y2 ) ( 5 - 2 ) 

with |#(0)) = * - It can be seen that V(-) G C a=1 (V(G) C U{G)) is linear 
and that cy := 2L 2 . Using Cea's lemma and a procedure similar to the followed 
in Chapter VII §5 in [5], it can be seen that the particular representation in (15. 2[) 
has approximation order v m = 2. Also from (|3.20[) we can observe that Hi := 
Ho+X 2 +Y 2 is self-adjoint, hence fflx is conservative and e ±l * Hl will be unitary with 
respecto to ||-|| H . If we take U to be given by the Crank- Nicholson approximation 
U := (1 + irHi)~ 1 (l — irM.%) wich is clearly a Pade approximant of e _4 * Hl and if 
we replace this approximation in (|4.38[) . since all the conditions are satisfied, then 
by C 14.21 and L 14.41 we will have that there exists an approximate solution 
to (|5.2|) that converges to the exact solution to (|5.1|) according to the estimate 

3 

\\\m)-\*i(t))\\ < c^ + ^h 3 T (5.3) 

where h G R + is the mesh size of the particular grid G\_h C G implemented. It is 
also important to note that if we denote the energy expectation evolution by {E) t 
and its discrete representation by (E}„ then it can be seen that 

ME>„ = ;((E)»+i-(E>J (5-4) 
= I((U" +1 * | Hi |U" +1 *o> - (U"*o| Hi |U n * » (5.5) 

T 

= (-(U-l)U n * |Hi|U n *o) 

T 

+ (U"* |H 1 |-(U-l)U n * > 

T 

+ ((U-l)U"*o|Hi |-(U- 1)U"* ) (5.6) 

T 

now it can be seen that the total variation for t = gives 

S (E)„ = i (U n * | [H l7 Hi] |U"* ) = (5.7) 
wich means that Sq kills (E) n , Vn, i.e., (E) is a discrete constant of motion. 
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Example 5.2. Nonlinear Schrddinger Equation. For quantum dynamical systems 
described by a nonlinear Schrddinger equation of the form 

i§j\m)--=*2\m)+*w)\ n w)) (5-8) 

with IV'(O)) := \ipo)) 2 < n G Z + and a G C. If we have that the quantum 
dynamical systems in this example evolve in the same media described in example I 
and have the same physical localization restrictions for its representative particles. 
Then the only condition that we need to check is the locally Lipschitz condition of 
V(X) := a\X\ n X, so we will need to restrict our attention to functions on a set 
S(G) C H that is caracterized by S{G) := {<j> G U{G) : \\<f>\\ < r, R+ 3 r < oo}, 
now, since \\X\ n X - \Y\ n Y\ < C„(|X|™ + \Y\ n )\X - Y\, where C n > is a value 
that depends on n, we will have that cy := 2\a\C n r n is the Lipschitz constant of 
V(-) with respect to W'Wy^/Qy If we use the same spatial particular projection than 
in the first example reewriting (|5.8p in the form 

i— := M" 1 P*M'P + a|*(i)|" (5.9) 

and if we also approximate U using the Crank-Nicholson approximation then by 
C \4-2\ and T \4-3\ we will have that there exists an approximate solution \i&i(t)) to 
problems of the form (|5.9j) that converges to the exact solution to ()5.1j) according 
to the estimate 

9\ty\f! r n+1 T 

WW)) - \Mt))\\c<X-r,rWC)) ^ ^ h + l-2| a "|C„r»r (2 + |a|r " r) (5 - 10) 

where h G R + is the mesh size of the particular grid Gij t G G implemented. 

In a similar way we can implement the operator techniques used here to study 
convergence and other related properties for other particular projectors with differ- 
ent approximation orders, also other Pade approximants than the Crank-Nicholson 
scheme can be used for the approximation of the integrating factor. 

6. Conclusion 

Particular projection methods can be implemented in the functional numerical 
analysis of several types of Schrddinger evolution equations, in this work we intro- 
duce some of the strategies that can be used and also we derive some estimates 
that can be very useful when we deal with the numerical solution of such important 
equations in quantum physics. 
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Appendix A. Some Results from Functional Analysis 

In this section we present some important results from functional analysis that 
are used in this work. 

Definition A.l. An operator T £ B(X) with X a Banach space, for wich 

\\T(x)-T(y)\\ < \\x-y\\,x,yeX (A.l) 

is called a contraction. If there is a K < 1 for wich \\T(x) — T(y)\\ < K \\x — y\\, T 
is called a strict contraction. 

Theorem A.l. Contraction Mapping Principle. A strict contraction T £ B(X) 
on a Banach space X has a unique fixed point, ie., there exists a unique x £ X 
such that T{x) = x. 

Definition A. 2. Succesive approximation methods for fixed points. Given a strict 
contraction T £ «Sf (X) on a banach space X , the sequence {x n } defined for some 
y £ X in the form 

Xn:= {T{xJl),n>l (A - 2) 
will be called an approximation sequence to x, and the iterative method will receive 
the name of successive approximations method. 

Remark A.l. A fixed point x £ X of a contraction T £ B(X) in a Banach space 
X , satisfies the following estimate 

\\x - x m \\ < K m (l - K)- 1 \\ Xl - x \\ . (A.3) 
with respect to its approximating sequence. 
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Definition A. 3. Locally Lipschitz Functions. A function f is locally Lipschitz 
,i.e., f € C a=1 (B r (0)), with B r (0) a closed ball of radius r centered in 0, if for 
< e r G R that depends on r, small enough, there exists c/ < oo with 

ll/(«) - f( v )\\ H ™(G) < c f \\ u ~ v\\ Hm{G) (A.4) 
w/ien 1 1 — v\\ Hm , G >. < e r . 

Definition A.4. Dissipative Operator. An elliptic operator A 6 Jtf(H n (G)) is said 
to be dissipative if we have that 

Re(Ax,x) < ,x e dom(A). (A. 5) 

in particular if equality holds, then A is said to be conservative. 

Theorem A. 2. If A £ £(H(G)) where H{G) is a discretizable Hilbert space, and 
if A is closed, densely defined and dissipative then it generates a contractive semi- 
group. 

Theorem A. 3. Stone's Theorem. A densely defined operator iA G Jz?(%) in a 
complex Hilbert space H is the generator of a strongly continuous unitary group on 
H if and only if A is self-adjoint. 
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